Nonlinear perturbations of the Kaluza-Klein monopole 
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We consider the nonlinear stability of the Kaluza-Klein monopole viewed as the static solution of 
the five dimensional vacuum Einstein equations. Using both numerical and analytical methods we 
give evidence that the Kaluza-Klein monopole is asymptotically stable within the cohomogeneity- 
two biaxial Bianchi IX ansatz recently introduced in Q- We also show that for sufficiently large 
perturbations the Kaluza-Klein monopole loses stability and collapses to a Kaluza-Klein black hole. 
The relevance of our results for the stability of BPS states in M/String theory is briefly discussed. 



INTRODUCTION 



The Taub-NUT instanton is the complete Riemannian 
Ricci-flat 4-manifold with the metric 
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where p > 0, < 6 < it, < (p < 2n, < %b < Air and 
m > is a free parameter which sets the scale. The Taub- 
NUT instanton is topologically R 4 with the surfaces of 
constant p being the squashed three-spheres. Although 
complete, the metric Q is not asymptotically Euclidean 
because the S 1 fibers approach a constant circumference 
at infinity (rather than growing like p). This kind of 
asymptotic behavior is called asymptotically locally flat. 
By adding the trivial term — dt 2 one obtains from (JIJ 
the regular static solution of the five dimensional vacuum 
Einstein equations 
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monopole falls within the BCS ansatz, the question of 
its nonlinear stability now becomes tractable and is the 
subject of this paper. 

The KK monopole plays an important role in M/String 
theory. For example, its metric product with six flat Eu- 
clidean dimensions gives a Ricci-flat eleven dimensional 
Lorentzian metric, which when reduced to ten spacetime 
dimensions may be interpreted as a D6-brane solution of 
Type HA String theory |4|. As such, it is an example of 
what is called a supersymmetric or BPS solution. BPS 
solutions of supergravity theories are widely believed by 
many string theorists to be absolutely stable, both classi- 
cally and quantum mechanically. However until now, this 
belief has never been tested in a fully non-linear setting. 
As we shall see, our results show that the classical stablity 
is not absolute since for sufficiently large perturbations a 
gravitational collapse to a black hole is possible. The im- 
plications of our results for BPS states will be discussed 
briefly in the conclusions. 

PRELIMINARIES 

) The BCS ansatz has the form |3] 
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This metric has been used in the past to obtain the four- 
dimensional magnetic monopole via a Kaluza-Klein re- 
duction along the coordinate ip |2| and since then it is 
usually called the Kaluza-Klein (KK) monopole. To the 
best of our knowledge, the role of the KK monopole in 
the dynamics of non-asymptotically flat initial data, and 
in particular its nonlinear stability, has not been stud- 
ied, probably because this problem appeared to require 
an analysis of the full five dimensional Einstein equations 
which is beyond currently available numerical tools, let 
alone the analytic ones. This situation has changed re- 
cently with the introduction of a new cohomogeneity-two 
symmetry reduction of the five dimensional vacuum Ein- 
stein equations (referred to below as the BCS ansatz |2j) 
which provided a framework for investigating the dynam- 
ics in a simple 1 + 1 dimensional setting. Since the KK 



where A, 8, and B are functions of t and r. In contrast 
to spherical symmetry, this ansatz possesses a dynami- 
cal degree of freedom, the field B, which measures the 
deformation of sphericity of the angular part of the met- 
ric. Substituting the metric J3J) into the vacuum Einstein 
equations one gets the following system of equations 
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where primes and dots denote derivatives with respect to 
r and t, respectively. The initial value problem for this 
system was investigated in 3] in the case of asymptot- 
ically flat initial data. Here we consider asymptotically 
locally flat solutions. 

The KK monopole (J3J) expressed in terms of the ansatz 
(|3j) takes the form 



m 



(1 



Ap 



B ° = I ^ f 1 + 3) . A) = ' 3 m '„, .<-"" = - 1, , . ( •". I 



(1 + ^) 8 / 3 
to 



where 



r» 1/2 1/2 /-, 1 P \ 
2m ' p ' I 1 H 

V TO/ 



1/6 



(6) 



Note that scaling invariance excludes the existence of 
nontrivial solutions which are everywhere regular and 
asymptotically flat. The KK monopole is the unique non- 
trivial regular static solution of the system (4). To see 
this note that static solutions which are analytic at the 
origin form a one-parameter family with the following 
asymptotic behavior for r — > 
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where b is the shooting parameter. Since the system (4) 
is scale invariant, it follows from the uniqueness of solu- 
tions of ordinary differential equations that, up to scaling, 
there is only one static solution which is analytic at the 
origin. It is easy to verify that the solution |5j) satisfies 
the regularity conditions Q with b — 1/(12to 2 ). For 
convenience, hereafter we choose units in which to = 1. 

LINEAR STABILITY 

Before discussing numerical results, we first demon- 
strate that the KK monopole is linearly stable within our 
ansatz. To this end, following the standard procedure we 
seek solutions of the system (4) in the form 

B(t,r) = B { P ) + B 1 {t,p), A{t,r)=A (p)+A 1 (t,p), 
6(t,r) = S (p) + S 1 (t,p), (8) 

where p is a function of r determined implicitly by l|^)l. 
Substituting the expansion JSJ into the system (4), lin- 
earizing and separating the time dependence B\(t,p) = 
cxp(— i\t)v\(p), after a straightforward but tedious cal- 
culation we get the eigenvalue equation for the spectrum 
of small perturbations around the KK monopole 
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Although this equation cannot be solved analytically in 
general, there is an explicit solution for A = 



vo 



P 



3 + 4p 



(10) 



This is the zero mode corresponding to the scaling free- 
dom - it can be obtained by the action of the scaling 
generator r-^ on B (p(r)). The zero mode l(TU)l has no ze- 
ros which implies by the standard Sturm-Liouville theory 
that the operator A has no negative eigenvalues. Thus, 
the KK monopole is linearly stable within our ansatz. 
Note that the zero mode is not a genuine eigenfunction 
because it is not square integrable. More precisely, on the 
Hilbert space L 2 ([0,oo),p(l + p)dp) the self-adjoint op- 
erator A has the purely continuous spectrum A 2 £ [0, 00). 
Nota bene, the Taub-NUT instanton Q is known to be 
linearly stable against transverse traceless perturbations. 
This follows from the work of Hawking and Pope J5| , who 
showed, using two linearly independent covariantly con- 
stant spinors which the Taub-NUT metric admits, that 
the spectrum of the Lichnerowicz operator, which gov- 
erns the linearized transverse traceless perturbations, is 
the same (apart from zero modes) as for the scalar Lapla- 
cian, and therefore non-negative. 

NUMERICAL RESULTS 

Using a fourth-order accurate finite difference code we 
have solved the system (4) numerically for several families 
of regular initial data which represent various perturba- 
tions of the KK monopole. The overall picture does not 
depend on the specific choice of a family. The results 
shown below were produced for initial data of the form 
(using the momentum variable P = e s A~ 1 B) 

B(0,r) = Bo(p(r)), P(0,r) =p g)'^, (11) 

where the amplitude p was varied and the parameters R 
and s were kept fixed. Note that although the perturba- 
tion Ulljl is exponentially localized, the induced pertur- 
bation of the function A has a 1/r 2 tail as follows from 
the hamiltonian constraint (4a). 

Small perturbations 

We have found that for small perturbations, that is for 
small values of the control parameter p, the perturba- 
tion is scattered off to infinity and the solution returns 
to equilibrium, i.e., it settles down to the KK monopole 
with the same parameter to. This is shown in Fig. 1. 
One might wonder why the parameter m does not change 
under perturbation. The reason is simple: the perturba- 
tions considered by us have finite energy while a change of 
to would require infinite energy. Here by energy we mean 
the energy measured with respect to a given KK refer- 
ence background as described by Deser and Soldate |||. 
This energy is determined by the next to leading order 
term in the asymptotic expansion at spatial infinity. 

The qualitative picture of pointwise convergence to the 
KK monopole is shown in Fig. 2. The details of this 
process, in particular the role of quasinormal modes and 
the rate of decay of a tail, will be pursued elsewhere. 

We have found that a sufficiently strong kick may 
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FIG. 1: Asymptotic stability of the Kaluza-Klein monopole. 
For initial data 1111 with a small amplitude (p = 0.1, R — 
3, s = 1) we plot a series of snapshots of the function B(t, r) 
(where t is central proper time). The dashed line shows the 
unperturbed KK monopole. During the evolution the excess 
energy of the perturbation is clearly seen to be radiated away 
to infinity and the solution returns to equilibrium. 



equations in five dimensions |?J 
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where 



A = p+3M-y / M 2 + 2P 2 , £ = p+M+y/M 2 + 2P 2 . 

Here p > 0, the parameter M is positive and the pa- 
rameter P (usually called the magnetic charge) has the 
range < P < 2M. For P < 2M the metric (O has the 
regular horizon at p — 0, so after |7j we shall refer to this 
solution as the Kaluza-Klein black hole. For P = 2M, 
the horizon degenerates to a point and one gets the KK 
monopole @ with m = 4M. 

Translating (12) into the ansatz 10 we obtain 
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FIG. 2: The convergence to the Kaluza-Klein monopole. For 
the same initial data as in Fig. 1 we plot the time series 
In |P(£,ro)| at ro = 2. The asymptotic power-law tail is seen 
at late times. 



Note that the leading order asymptotic behavior for large 
r does not depend on M and is determined only by P 
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Although the KK black holes are known in closed form, 
for our purposes it is helpful to look at them from the 
viewpoint of the shooting procedure. Assuming that 
there is a non-degenerate horizon at r — rjj, we obtain 
from (4) the following behavior 
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hence for each rjj > there is a one-parameter family 
of local solutions with a regular horizon parametrized by 
(3 = B(rH) > 0. Shooting these local solutions towards 
infinity we get the KK black holes l|13|) where 



destabilize the KK monopole. In order to understand 
what is the fate of a strongly perturbed KK monopole 
we need to make a digression. 

Kaluza-Klein black holes 

The KK monopole J2J is in fact a special case of the 
two-parameter family of solutions of the vacuum Einstein 
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Large perturbations 

Now, we return to the discussion of numerical results. 
We have found that for sufficiently large perturbations 
the inner part of the KK monopole collapses and a hori- 
zon forms at some r — rn > 0. Outside the horizon the 
solution settles down to the KK black hole (see Fig. 3). 



As we pointed out above, the perturbations considered 
by us have finite energy and therefore they do not change 
the leading order asymptotic behavior (fl 51) . In this sense 
the magnetic charge P can be viewed as the constant of 
motion. For the KK monopole P = 2M = m/2, so all our 
configurations have the same magnetic charge P = 1/2 
(recall that we use units in which m = 1) and the end- 
states of evolution differ only by M. 
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FIG. 3: Instability of the Kaluza-Klein monopole for large 
perturbations. For initial data (lilt with a large amplitude 
(p — 0.3, R — 3, s = 1) we plot a series of snapshots of the 
function A(t,r) (where t is central proper time). During the 
evolution A(t, r) drops to zero at r — th ~ 1-68 which signals 
the formation of an apparent horizon there. Outside the hori- 
zon the solution relaxes to the Kaluza-Klein black hole (I13H 
(shown by the dashed line) with M — 0.502. 

As we approach the threshold of collapse from above, 
the parameter M decreases towards the limiting value 
M = 1/4 corresponding to the KK monopole. At the 
same time, as follows from (|17f) . both j3 and r# tend 
to zero. In this limit, close to the horizon the KK black 
hole is very well approximated by the Schwarzschild black 
hole. Near the threshold we observe the type II discretely 
self-similar critical behavior governed by the same critical 
solution as in the asymptotically flat situation [3j (see 
Fig. 4). This is not surprising in view of the fact that 
critical collapse is a local phenomenon, and, as such, does 
not depend on the far field behavior. 

CONCLUSIONS 

We have shown that although it is classically stable 
against small perturbations, the KK monopole is clas- 
sically unstable against perturbations which are large 
enough to allow gravitational collapse to form a black 
hole. Classically, the black hole is probably stable. Quan- 
tum mechanically one expects it to emit Hawking radia- 
tion which can carry off energy but not magnetic charge. 
In the long run therefore one might expect the system 
to settle down to the original KK monopole. If this is 
so, then the common intuition about BPS states would 



be correct, provided one bears in mind that the stability 
can only hold by virtue of quantum mechanical effects. 
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FIG. 4: Critical behavior. For supercritical solutions we plot 
the horizon radius vs. the distance to the threshold on the log- 
log scale. The fit of the power law ln(r_tf) = (7/2) ln(p — p*) + 
const yields 7 w 0.33 which agrees (up to numerical errors) 
with the critical exponent obtained in Q- The echoing period 
A S3 0.47 (read off from the period of wiggles around the linear 
fit above and computed independently from the spatial shift 
between the echoes) is also the same. 

Finally, we remark that, besides the KK monopole, the 
BCS ansatz incorporates other non-asymptotically flat 
solutions. In particular, an exact regular time dependent 
solution is known [8j 
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This solution can be viewed as the time evolving KK 
monopole whose S 1 fibers pinch off to a point for t — > 
00. It would be interesting to determine the role of this 
solution in dynamics. 
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